Three-body Halos. 
IV. Momentum Distributions after Fragmentation 



E. Garrido, D.V. Fedorov and A.S. Jensen 
Institute of Physics and Astronomy, 
Aarhus University, DK-8000 Aarhus C, Denmark 
(February 9, 2008) 



I. INTRODUCTION 



Momentum distributions of particles from nuclear break-up of fast three-body halos are calculated 
incorporating effects of final state interactions. The same two-body interactions between the particles 
\& ■ are used to calculate both the ground state structure and the final state of the reaction processes. The 

0^ ' ground state wave function reproduces size and energy of the halo nucleus. First we give a general 

and detailed description of the method. Then we investigate the effect of final state interactions in 
neutron removal and core break-up reactions for one- and two-dimensional momentum distributions. 
(3JT)' We compute specifically core and neutron momentum distributions from Li fragmentation and 

compare those with available experimental data. We conclude that 11 Li must have a neutron-core 
relative state containing a p-state admixture of 20%-30%. The fragmentation data also strongly 
suggest that 10 Li has an s-state at about 50 keV and a p-state around 500 keV. 
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In the latest years a big effort has been made to investigate one of the most striking new features in nuclear physics: 
halo nuclei. They are weakly bound and spatially extended systems where one or two particles have a high probability 
of being at distances larger than the typical nuclear radius. General properties and the criteria for the ocurrence of 
■ nuclear halos have been discussed in [jJ-Q] and a review of nuclear halos has recently been published 
"^T* | The appearence of the halo structure in light nuclei along the neutron dripline Q gives special interest to this 
region of nuclei. The most general properties of them are well described by few-body models, dividing the degrees of 
freedom into the approximately frozen (core) and the active (halo) degrees of freedom. In particular, special attention 
has been paid to Borromean systems consisting of three-body bound systems, where all the two-body subsystems are 
unbound |]||. The prototypes nuclear halos are 6 He ( 4 He+n+n) and n Li ( 9 Li+n+n), both thoroughly discussed in 
• t-h , a general theoretical framework in Q . 

The ability to produce secondary beams of halo nuclei opens the possibility of investigating their structure by 
measuring the momentum distribution of "particles" resulting from fragmentation reactions pToj— pTq] . However one 
major problem in the interpretation of such measurements is the inherent mixture of effects from the original structure 
of the projectile and the reaction mechanism. Both should therefore be properly incorporated in model calculations. 

The simplest picture used to describe this kind of fragmentation reactions is formulated in the sudden approximation, 
where one of the three particles in the projectile is instantaneously removed by the target, while the other two particles 
remain undisturbed. Clearly this can only be justified for reaction times much shorter than the characteristic time for 
the motion of the three particles in the system. Since the projectile is weakly bound this requirement is well fulfilled 
for any high-energy beam. The observed momentum distribution should then provide direct information about the 
three-body wave function. The validity of this model as a first description of the reaction has been discussed in the 
litterature, see for example p7 - 19 . 



The connection between the wave function and the momentum distributions is more complicated than originally 
anticipated when the early experiments were designed. It was found recently that momentum distributions are highly 
affected by the final state interaction (FSI) between the two particles remaining after the collision, especially when 
low- lying resonances are present |TJ,|l8|,^0| . This essential ingredient must then be incorporated into the model. The 
final state two-body interaction is at the same time determining the three-body structure of the projectile (halo 
nucleus), and therefore a consistent treatment of the initial and final states is needed to allow reliable interpretations 
of the experimental data. 

This paper is number four in a series discussing the general properties of three-body halo systems. In the first 
of these || we discussed Borromean systems, in the second [|J we extended the discussion to general three-body 
systems and in the third we discussed the effects of finite spin of the particles. In this paper we investigate 
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fragmentation reactions of Borromean projectiles by use of the sudden approximation. A brief report containing a 
few of the pertinent results for 11 Li is published p2| . Our main purpose is here to formulate a general procedure to 
incorporate the effects of final state two-body interactions. This two-body interaction is consistently included both 
in the final state wave function and in the description of the three-body projectile. We apply the method in detailed 
analyses of measured n Li fragmentation reactions. We thereby obtain essential information about the unbound 10 Li 
system. 

The paper is organized as follows: The mathematical formalism is developed in Section 2, which to facilitate the 
understanding is divided in five subsections. The general numerical results are presented in Section 3, where we as 
illustration use parameters relevant for n Li fragmentation reactions. In Section 4 we perform realistic calculations 
for 11 Li and compare to experimental data. Finally, Section 5 contains a summary and the conclusions. Some of the 
mathematical definitions are collected in an appendix. 

II. METHOD 

To describe three-body halo fragmentation reactions we shall use the sudden approximation. We are then assuming 
a process where a high-energy three-body halo projectile instantaneously looses one of the particles without disturbing 
the remaining two. We also assume a light target and we shall therefore not consider Coulomb dissociation process, 
which then only contribute marginally . 

We work in the center of mass system of the three-body projectile and denote by k and q the total and relative 
momentum of the two remaining particles in the final state. The transition matrix of the reaction in the sudden 
approximation is then given by 

M(k,q) oc (e ik ' R e iq - r |*> , (1) 

where VP is the three-body wave function, r is the distance between the two reamining particles, and R is the distance 
between the center of mass of the two-body system and the removed particle. 

This transition matrix describes the final state as a plane wave, which means that the interactions between the 
particles are neglected. In principle all of them should be included in the calculation, but under the experimental 
setup where one of the particles suddenly is removed, only the interaction between the two non-disturbed particles 
is present. To account for this final state two-body interaction the plane wave e lq r has to be replaced by the 
corresponding distorted two-body wave function u>(q, r). 

In the following subsections we describe how we construct the initial three-body halo wave function \& as well as the 
distorted two-body wave function w(q, r). Finally we compute the transtion matrix and the momentum distributions. 

A. Initial Three-body Wave Function 

We shall use the Jacobi coordinates, basically defined as the relative coordinates between two of the particles (x) 
and between their center of mass and the third particle (y). Three sets of Jacobi coordinates can be constructed 
and the precise definitions and the corresponding sets of hyperspherical coordinates (p, a, Cl x , fl y ) are given in the 
appendix. One of the three sets of hyperspherical coordinates is sufficient for a complete description of the system, 
but the subtle correlations are easier to describe by use of all three sets. The volume element is given by p 5 dfldp, 
where d£l — sin 2 a cos 2 adadSl x dSl v . 

The total wave function jm of the three-body system (with total spin J and projection M) is written as a sum 
of three components tpjM > which m turn for each p are expanded in a complete set of generalized angular functions 

3 3 

*«=E^ji) = -^E X>Sm(/>^) . (2) 

i— 1 ' n i— 1 

where the radial expansion coefficients f n (p) are component independent and p -5 / 2 is the phase space factor. 
These wave functions satisfy the three Faddeev equations p^j24] | 

(T - E)i>f M + V jk {i,f M + i,f M + ^f M ) = , (3) 

where Vjk are the two-body interactions, E is the total energy, T is the kinetic energy operator and {i, j, k} is a cyclic 
permutation of {1,2,3}. 
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The angular functions are now chosen for each p as the eigenfunctions of the angular part of the Faddeev equations 



where A:} again is a cyclic permutation of {1,2,3}, m is an arbitrary normalization mass, and A 2 is the p- 
independent part of the kinetic energy operator. The analytic expressions for A 2 and the kinetic energy operator can 
for instance be found in 0. 

The angular functions ^^j M (p,rti) are expanded in terms of the complete set of hyperspherical harmonics 
Yfg (ofj , Cl Xi , il Vi ) , where the quantum number K usually is called the hypermoment, £ x and t y are the orbital 
angular momenta associated with x and y, and L is the coupling of these angular momenta. The result is 

€jmM) = E ^t^r.. m sO>)[rtXte>to*M®x^ a ] JM , (5) 

(i) 

where C are expansion coefhcients and x s s s ^ s ^ ne three-body spin function. The spins of the two particles, connected 
by the x coordinate, couple to the spin s x , which coupled to the spin s y of the third particle results in the total spin 
S of the three-body system. 

The radial expansion coefficients f n (p) are obtained from a coupled set of "radial" differential equations B, i.e. 



d 2 2mE 



\ ( A„(P) + /"(P) + E {- 2P nn>^- - QnrA fn'(p) = , (6) 



dp 2 ti 2 p 2 V w ' 4JJ dp 

TV 

where the functions P and Q are defined as angular integrals: 

3 

P nn ,{p) = E / dQ^*(p,Q)^)(p,n), 

i,j = l 

Q nn ,( P ) = E / dn*w*( P ,n)^^( P) n) . 

Finally from eqs. (|^) and (|^) we obtain the expression of the three-body wave function 



(7) 



(8) 



JM 



*^=-^72EU°)E E ^^ y ^s(p)[Y^(a 4 ,^,^)® x ^ s j ■ , (9) 
where the radial function f n (p) and the coefficients C^ Ki t Ls s (p) must be computed numerically. 

B. Final-state Two-body Wave Function 

In the fragmentation reaction the target suddenly removes one of the particles from the projectile and the remaining 
interaction in the final state therefore acts between the two non-disturbed particles. Assuming that the two-body 
interaction does not mix two-body states with different spin s x and relative orbital angular momentum £ x , we can 
expand the two-body wave function in partial waves, see p5J, as 



^(k X)X ) = yjll- J2 utjk x ,x)yzi sx (n x ) 



x 22 (Zxrn,£ x s x a x \j x m x )i tx Y ia:m ^(Q, k3i ) , (10) 

mi a> =— ia, 

where u x is the spin projection of s x , u is the radial and y the angular distorted wave function. The angles fl x and 
f2fc x define the direction of x and k. x . These assumptions are usually strictly valid. The only exception arises from 
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the tensor interaction and even then the resulting mixing is often very small and therefore insignificant in the present 
context. 

The two-body space and momentum coordinates x and are defined consistently with the Jacobi coordinates 
introduced in the three-body wave function (see appendix) 




where r is the relative spatial coordinate and p r the relative momentum. The mass m is the arbitrary normalization 
mass and /i is the reduced mass of the two-body system. 

The radial functions u J f s (k x , x) are obtained numerically by solving the Schrodinger equation with the appropriate 

two-body potential V{x) 

|U , ,:/,,..,:, + (kl - |?Vt» - u£ Sx (k x ,x) = , (12) 

where 

vtsM = J '/<^j;1^j<t)r;xij-;;' ( ...io, i . ( i 3 ) 

When the interaction between the two particles is neglected the solution of eq.(|l^) is (k x x)j^ sl .{k x x), and the 
expansion in eq.(p^|) reduces to the usual expansion of a plane wave in terms of spherical Bessel functions ji. 



C. Transition Matrix 



To compute the transition matrix in eq. (|l]) it is convenient to rewrite the three-body wave function in eq. (|^) in terms 
of the set of Jacobi coordinates where x is related to the two particles remaining after the fragmentation. This means 
that all three Faddeev components in eq.(||) have to be expressed in terms of the chosen set of Jacobi coordinates. 
The three-body wave function can then be written 



#.™*(x,y) = -572 53^00 c nKlxl , vIjSx s{p)WfX{(x,&x,siy)®x 



1 JM 



(_s x s v S 



(14) 



where CnKt^tyLs^sip) are the coefficients of the expansion after the transformation. 

The calculation of the transition matrix involves the overlap of eq.(fL4|) and the function e* kH ' y w SxCTx (k K , x) (see 
cqs.(0) and ( |Io| ) and note that k • R = k y • y, where k y is the momentum related to y). After integrating analytically 
over £l x and Q y we get for each X n the following expression for the transition matrix: 



M 



JM 



^_ l yJ+2S-2M+l y +s y -s :I -t3 S 2'2~-2 



h 



JxJy y frl —m x — 1Tl v 



Ly Oy 



-m y rri£ y a y 



jx $x 

-m x a x 



J jx jy 

L £ X £y 
S S X Sy 



where a = \/2a + 1, a y is the projection of s yi ( ) and { } are the usual 3 J and 9 J symbols, see 



(15) 



We have 



introduced the hyperspherical coordinates in momentum space n = w/cjj; + ky, a K — arctan (k x /k y ) and the function 



7 l B f,,- ( K > a «) is g iven b y 



^x^xj 



/■CO 

^ ««) = ^ +IV E N K tV \ P b/2 dpf(p)C nK LA y Ls,s{p) 
K J ° 



(16) 



tt/2 



da{sin af*+ 2 {cos a) l y +2 Pt* +h '" y+h (cos(2a)) 3iy (k y y)^-uf s {k x ,x) 
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where v = (K — £ x — ly)/2, PJ° 2 ' v 2 is a Jacobi polynomial and 



JV K 



v\{ V + t x +l y + l)\2(K + 2) V /2 

i> + 4 + f)i> + ^ + §) J 



When the interaction between the two particles in the final state is ignored the integral over a in eq.([L6|) can be 
carried out analytically, i.e. 

7T/2 j j 

da(sina)^+ 2 (cosa)^+ 2 P^ + ^ + 2 (cos(2a)) j, y (^y)^ x (fei) = 



o 



x(-l)^(sina K )^(cosa K )^P,^ + ^ + Mcos(2a K ))^±^ (18) 



D. Momentum Distributions 

The cross section or momentum distribution is now obtained by squaring the transition matrix and subsequently 
averaging over initial states and summing over final states: 

«E E \K^ y (^,K)\ 2 - d9) 



The volume element dk x dk y is written 

d\s. x dk y = k x dk x L dk x difik :c kydkydrtk y , (20) 

where k x = fe x sin k\ = k x cos9h x , and ((9^ , ) are the polar and azimuthal angles of the vector k x . The 
differential cross section in e q.(|19| ) should be integrated over all unobserved variables. 

Substituting eq.(|T5|) into (|19|) we obtain, after analytical integration over fifc s and (fk x , the following expression for 
the three-dimensional differential cross section or momentum distribution: 

d3(J '„2 u± 2 V- V- V- / 2 ■ 2 t f, 6 S, rfyLS 



K ^xS E E E^ ^ LL'§§'I£ am .jK,a K )I l l amjm (K,a K ) 




dkydk^jr dkx ^ 

[ J jx jy I [ J jx jy I 

(21) 

By integrating numerically over k y and k x , we get the one-dimensional relative momentum (fcl) distribution of the 

remaining particles. By integrating over k y and k x , we get instead the two-dimensional relative momentum (k x ) 
distribution. 

It should be noted that we have not specified any coordinate system and the axes x, y, and z are therefore completely 
arbitrary. Thus, in the sudden approximation the longitudinal and transverse momentum distributions are identical. 

E. Transformation to the Center of Mass of the Three-body System 

Comparison with the measured momentum distributions requires transformation to the center of mass system of 
the projectile. To do this we construct the momentum p of one of the particles in the final state relative to the center 
of mass of the projectile as a linear combination of k x and k y 

p = dik x + b t k y . (22) 



If the Jacobi coordinate x refers to particles 1 and 2, then a and b take the set of values (see eqs.(A4) and ( A5)) 



. 1 mim 2 \ 1/2 l mi f 1 (mi + m 2 )m 3 \ 1/2 
at = ; , 01 = ; ; ; , (23) 



m mi + m 2 / m\ + m 2 \ m m\ + m 2 + tus 
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v m mi + 1^2 J mi + m 2 \ m m\ + m,2 + 

when we compute the momentum distributions of particle 1 and 2, respectively. 

Using the momentum p of one of the particles as the variable instead of we obtain the relation 

d 6 a 1 d 6 a 



oc 



E \K A L y , y (^K)\ 2 , (25) 



where a 3 arises from the Jacobi determinant for the transformation. 

As before we must integrate eq.(^5|) over the unobserved quantities, i.e. k y and some of the components of p. It is 
then convenient to use both the initial (arbitrary) coordinate system and a rotated system where the z-axis is along 
p. The volume element is then written as 

dpdky = p~ L dp ± dp^d(p p kydkyd£l' k ^ , (26) 

where p 1 - — ps\n8 p , pH = pcos# p , (9 p ,(p p ) define the direction of p in the initial coordinate system and the angles 
Q' k — (0' k , ip' k ) give the direction of k y in the rotated system. 

We express the transition matrix in eq. (|l5| ) as a function of p referred to the initial system and k y referred to the 
rotated system. Then we need the relations 

Y em (9 ky ,vk v ) =Y,vi m ,(v P ,0 P ,0)Y em <(9' kv ^' ky ) , (27) 

m' 

Y im {6 kx , Vkx ) = J2Km'(v P ,0 P ,0)Y em ,(9l^' k J , (28) 

where T>(ip p ,9 p ,0) is the rotation matrix (the P-functions) defined in |2q] , and the angles [9' k ,ip' k ) defining the 
direction of k^ in the rotated system are given by 

p — bik v cos Q'u 

cos 9' kx = \ * , (29) 

Q>i Kx 

, _ { <p' ky if a i/h < ,» n x 

^~\^ h +7T if CLi/biX) [M) 

Finally, we need to express k and a K in terms of the new variables p, k y and 9' k , i.e. 
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« 2 = ^2 [P 2 + b 1 k l ~ 2h ^ k V C0S d 'k v ) + K > ( 31 ) 

kx (p^ + blkl-2hpk y cos9' ky )^ n 
tana K = — = — * , < a K < - . (32) 

The numerically obtained function I t y ^ . (n,a R ) in eq. (p5|) depends on 9' k and only the integrations over ip p and 
tp' k in eq. ( p5| ) can be done analytically. We then obtain the four-dimensional differential cross section or momentum 
distribution relative to the center of mass of the projectile 

— f* ociyi- s yi y y y (_!)•.— 

dkydO', dp ± dpW yF afn 2 "y ^ ' 

^SxJxLS e' m j' m L<S< lyjyl'j' 



J Jx jy 

S Sx Sy 



{ji \ f 2 { $x ^ Lx \ ( £y t-'y L x \ J j y j' y L x 1 

s> si i/tr o o ; ^ o o o ){e>iy s y \ 



G 



where we assumed that <ij < 0. When a,; > 0, the summation over N x should include an extra factor (— 1) N ". 

Three integrations must be done numerically in order to get the one-dimensional (pH) and two-dimensional (jp^) 
momentum distributions of one of the particles in the final state relative to the center of mass of the projectile. 

If we assume that one of the particles in the final state, say particle 1, has infinite mass, then the momentum 
distribution of particle 2 relative to the three-body center of mass should coincide with the momentum distribution 
relative to particle 1. In other words, eq.(pl|) should be recovered by integrating eq.(|33|) over 6' k . Then 6, — (see 



([3 2D) and the integration over 6' k can easily be done analytically leading directly to eq.(|2 



cq.(|2J)), p = aik x and therefore 0' k = and P L X (1) = 6if x ,o- Now k and a K are independent of 9' k ^ (eqs.(31) and 



III. NUMERICAL RESULTS 



In this section we apply the method to the fragmentation of two-neutron halo nuclei (core+n+n) . Two different 
kinds of processes will be considered: neutron removal and core break-up reactions. In the first case one of the 
neutrons is removed by the target, and the neutron-core interaction in the final state must be considered. In the 
second case, where the core is violently removed from the projectile, the neutron-neutron interaction is involved in 
the final state. 

We take parameters corresponding to n Li in our examples. We first specify the two-body potentials used in the 
calculation. We then investigate the momentum distributions and how they are affected by final state interactions 
and the structure of the projectile. 



A. Two-body Potentials 



Loosely bound systems are mainly sensitive to the low-energy properties of the potentials. We shall therefore use 
relatively simple potentials reproducing the available low-energy scattering data and still allowing extensive three-body 
calculations. 

For the neutron-neutron interaction in the singlet s-wave we use the gaussian potential introduced in ref. Q . The 
extension to the triplet p-wave is made by including spin-spin, spin-orbit, and tensor terms. The neutron-neutron 
potential has then the form 

V nn = (V c + V ss s nl ■ s„ 2 + VrS±2 + V so \ nn ■ s nn J exp [-(r/&„„) 2 ] , (34) 

where s„i and s„2 are the spins of the two neutrons, s„„ = s„i + s„2, l nn is the relative neutron-neutron orbital 
angular momentum and S\2 is the ususal tensor operator. 

The strength parameters V c , V ss , Vt, and V so , and the range parameter b nn are adjusted to reproduce the following 
scattering lengths a and the s-wave effective range r e |27j] : 

aCS ) = 18.8 fin ,r e ( x S ) = 2.76 fin , (35) 
a( 3 P ) = 3.6 fin , a( 3 Pi) = -2.0 fin , a( 3 P 2 ) = 0.30 fin . (36) 

We then obtain values for the parameters of the potential 

V c = 2.92 MeV , V ss = 45.22 MeV , (37) 
V T = 26.85 MeV , V so = -12.08 MeV , (38) 
6„ n = 1.8ftn . (39) 

For the neutron-core potential we also assume a gaussian shape and we use the parametrization 

Vlc ] =V S (1 + j s Sc ■ s„) exp [-{r/b nc f\ , (40) 
V^J = (V t + VjXc ■ s«c) exp [~{r/b nc f] , (41) 

where s c and s„ are the spins of the core and the neutron, s nc = s c + s„ and \ nc is the relative neutron-core orbital 
angular momentum. This choice is convenient when we adjust the all-decisive energy positions of the virtual states 
and the resonances in the neutron-core subsystem. The spin splitting term s c • s„ is for simplicity only introduced in 
the s-state, but it could as well be included for other partial waves. 
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FIG. 1. One-dimensional core momentum distribution for a neutron removal process in n Li fragmentation. Thin curves 
correspond to calculations without FSI, while FSI are included for the thick curves. In the upper part the virtual s^-state 
is placed at 100 keV (solid curve), 200 keV (short-dashed curve), and 300 keV (long-dashed curve), respectively. In the lower 
part the p!/ 2 -resonance is placed at 100 keV (solid curve), 200 keV (short-dashed curve) and 300 keV (long-dashed curve), 
respectively. The momentum of the core is referred to the center of mass of the Li projectile. The Li spin is assumed to be 
zero. 
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TABLE I. Energy positions of the virtual s 1(/2 -state and the p 1(/2 -resonance in the neutron- 9 Li subsystem (first and second 
column) and the parameters of the neutron-core potential for each case. The last two columns give the s- and p-wave content 
in the neutron- 9 Li subsystem in the three-body wave function of n Li. 
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B. Momentum Distributions for Spin Zero Core 



We analyze the general properties of momentum distributions and the effects of final state interactions by taking 
11 Li as an example. We include s- and p- waves in the Faddeev components and neglect all higher orbital angular 
momenta, since their contributions are expected to be exceedingly small. 

We first consider the simpler case where the spin of the 9 Li-core is assumed to be zero. The spin quantum numbers 
of 10 Li and n Li are correspondingly 1/2 and 0. The range parameter b nc in eqs. ( ff0|) and (El]) is assumed to be 2.55 
fm ||. The strength parameter V s is used to modify the position of the virtual Si/ 2 -state in "*Li. The parameter Vi=± 
is chosen as —4.5 MeV and the spin-orbit parameter V so must be large to reproduce the experimental binding energy 
(295±35 keV ^8|) and root mean square radius (3.1±0.3 fm of 11 Li. With this choice of the parameters we have 
a simple structure for the three-body system, where only the positions of the virtual Si/2-state and the p^-resonance 
in 10 Li are relevant. In particular, then a p 3 / 2 -resonance appears at a very large energy. We can now investigate 
how the momentum distributions are influenced by final state interactions for different structures of the neutron-core 
subsystem. 

Comparing the cases of low-lying virtual s-states and low-lying p-resonances are of special interest. Six different 
situations are considered: The virtual Si/2-state in 10 Li is at 100, 200, and 300 keV, and the pi/2-resonance is at 100, 
200, and 300 keV. Because of the requirement of reproducing the experimental binding energy and rms radius in n Li 
both the virtual s^-state and the px/2-resonance can not simultaneously be at a low energy In table I we show 
the values of the parameters used in the neutron-core interaction as well as the resulting s- and p-wave content in the 
neutron-core subsystem of the total three-body wave function. 

In fig.|l| we show the one-dimensional core momentum distribution for a neutron removal reaction measured in the 
center of mass of the three-body projectile. To make the comparison easier all the curves have been scaled to the same 
maximum. The thin curves are the calculations without final state interactions while the thick curves are obtained 
with final state interactions. In the upper part of the figure the virtual sm-state has been placed at 100 keV (solid 
curve), 200 keV (short-dashed curve) and 300 keV (long-dashed curve), and in the lower part the p^-resonance 
energy is at 100 keV (solid curve), 200 keV (short-dashed curve) and 300 keV (long-dashed curve). The full width at 
half maximum (FWHM) for the momentum distributions in the figure are given in the second and third columns of 
table II. 

From fig.|l] and table II we observe the two important features. First, for a given neutron- 9 Li structure final state 
interactions always make the momentum distributions narrower. The change is not drastic, due to the large mass 
of the core, but it is still important. The largest effects appear for low-lying virtual Si/2-states, where the FWHM 
decreases by 15-19%. For low-lying p-resonances the FWHM decreases by less than 10%. Thus, the larger the s-wave 
content in the neutron-core subsystem, the larger the effect of final state interactions. 

Secondly, the width of the momentum distributions increases from 55 MeV/c, when the virtual s^-state is at 100 
keV, to 87 MeV/c, when the pi/2-resonance is at 100 keV. The corresponding numbers without final state interactions 
are 68 MeV/c and 92 MeV/c. Comparison with the experimental value in table II strongly suggests a low-lying 
virtual s-state in 10 Li. The corresponding s-wave content decreases from 97% for the narrowest to only 22% for the 
broadest distribution. Thus the larger the s-wave content in the initial neutron-core subsystem of 11 Li, the narrower 
the momentum distribution. 



Energy 
(keV) 




E, = 200 



Experimental 



TABLE II. Full width at half maximum in MeV/c of the one-dimensional core and neutron momentum distributions from 
Li fragmentation for different positions of the virtual Si/2-state and the p^-resonance. Columns 2 and 3 refer to core 
momentum distributions, while 4 to 6 refer to neutron momentum distributions. The label (n-c) indicates that final state 
interactions between 9 Li and the neutron have been included in this neutron removal process. The label (n-n) indicates that 
final state interactions between the two neutrons have been included (core break-up) [a] Data from pq |. [b] See, for instance, 
@. [c] Data from [jTJ. 



These two facts can be explained by the effects of the centrifugal barrier. A high p-wave content in the neutron-core 
subsystem increases the importance of the centrifugal barrier, see eq.((l^). Consequently the effects of distortion from 
the short-range potentials is smaller for p-waves than for s-waves. Therefore the effects of final state interactions 
decreases when the p-wave content increases. Furthermore, the probability for the neutron to be outside the range 
of the two-body potential drastically decreases for p- waves This decreases the spatial extension of the initial 

system when the p-wave content increases, and the momentum distribution is then broader. 




25 50 75 

p" (MeV/c) 

FIG. 2. The same as fig.[j] for one-dimensional neutron momentum distributions. 

In fig]^ we show the one-dimensional neutron momentum distributions for a neutron removal process. The notation 
to distinguish the curves is the same as in figj^. The FWHM for these distributions are shown in the fourth and fifth 
columns of tabic II. Let us focus first on the momentum distributions without final state interactions (thin curves). As 
in fig.0 we observe that a higher s-wave content in the neutron- 9 Li subsystem decreases the widths of the momentum 
distributions. The FWHM decreases from 94 MeV/c to 51 MeV/c when the s-wave content increases from 22% to 
97%. 

When final state interactions arc included in the calculations (thick curves), the change in the neutron momentum 
distributions is much more significant than for core momentum distributions. For low-lying virtual s^-states (upper 
part of fig.|^) the width of the distribution decreases for a virtual state at 100 keV by as much as 45% due to final 
state interactions. We still observe that the lower the virtual Si/2-state (high s-wave content in the neutron-core 
subsystem) the narrower the momentum distribution. For low- lying p^-resonances in 10 Li (lower part of fig.|^) we 
also observe a drastic modification of the neutron momentum distributions due to final state interactions. However, 
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we now see that the lower the p^-resonance (higher p-wave content) the narrower the momentum distribution. This 
result seems to contradict our experience from the core momentum distribution where the presence of the centrifugal 
barrier was decisive. But the effect of final state interactions on the much lighter neutron in a low-lying p-resonance 
is apparently much larger than the simultaneous broadening due to the centrifugal barrier. 

We can try to estimate the peak position of the p-wave contribution to the momentum distribution. The relative 

momentum (p r ) between the neutron and the core in the final state preferably takes the value *J 2fj,E TCS /h 2 , where 
E Tes is the energy of the resonance and /x is the reduced mass of the two-body system. The momentum of the neutron 
relative to the center of mass of the projectile is very close to k x , since the constant bi in eq.(|2^) is close to zero due 
to the large mass of the core. Then the neutron momentum distributions in fig.^j would be concentrated around the 
relative momentum p r = yf /i/mk x , see eq. (|Tl|) ■ If the resonance energy is low, the final state interactions would try 
to create a narrow momentum distribution. This argument is valid for both low-lying virtual Si/2-states and low- lying 
Pj/2-resonances. Thus, we can obtain a narrow neutron momentum distribution with FWHM compatible with the 
experimental value in table II both with low- lying S\/2 and pi/2-states. 

For two-dimensional neutron and core momentum distributions in a neutron removal process we obtain conclusions 
similar to those extracted from fig s.0 and § As a general rule we can say that momentum distributions are narrower 
and more affected by final state interactions when a low-lying virtual s-state is present in the neutron-core subsystem 
(high s-wave content) than with a low- lying p-resonance (low s-wave content). As for the one-dimensional distributions 
the neutron momentum distributions influenced by the final state interactions do not follow this rule. We can also 
get a narrow momentum distribution with a low- lying p^-rcsonancc in 10 Li. 
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FIG. 3. Two-dimensional neutron momentum distribution for a neutron removal process in Li fragmentation. FSI has been 
included in all the calculations. Solid curves represent the cases of low-lying virtual Si/2-states at 100 keV (thick curve) and 
200 keV (thin curve) in the neutron- 9 Li subsystem. Dashed curves represent the case of low-lying p^-resonances at 100 keV 
(thick curve) and 200 keV (thin curve) in the neutron- 9 Li subsystem. The momentum of the neutron is referred to the center 
of mass of the n Li projectile. The x Li spin is assumed to be zero. 

In fig|^ we show the two-dimensional neutron momentum distributions for a neutron removal process. The solid 
curves correspond to the cases where a virtual Si/ 2 -state is present at 100 keV (thick solid curve) and at 200 keV 
(thin solid curve) in the neutron-core subsystem, while the dashed curves represent the cases where a pi/2-resonance 
is at 100 keV (thick dashed curve) and at 200 keV (thin dashed curve). In all these cases final state interactions have 
been included in the calculation. 

For a neutron-core subsystem with a p-resonance at a given energy E ICS the momentum distribution (with FSI) 
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will exhibit a peak at the relative neutron-core momentum p r = y/ fi/mk x = y 2(iE ICS /h 2 . For E rea = 100 keV and 
200 keV, we then expect the peaks at k x ~ 14 MeV/c and at 19 MeV/c. In the two-dimensional neutron momentum 
distributions, we integrate away the dependence on k\ and refer to the center of mass of the projectile. Furthermore, 
even after including the s-wave contribution, which is 22% for the thick dashed curve and 36% for the thin dashed 
curve in fig.|], we still can see a clear bump in the momentum distributions corresponding to low-lying p-resonances. 
Thus, the effects of final state interactions on the two-dimensional neutron momentum distributions results in shapes, 
which depend very much on the structure of the low-lying state. 
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FIG. 4. The same as in figjl] for one-dimensional neutron momentum distributions in a core break-up reaction. 

Let us now turn to core break-up reactions, where the core is violently destroyed during the interaction, but where 
the sudden approximation still is valid. Then the final state interaction (neutron-neutron interaction) is well known 
and the uncertainties involved in the calculations are reduced. In fig^J we show one-dimensional neutron momentum 
distributions for core break-up reactions. The FWHM for the curves in the figure is given in the fourth and sixth 
columns of table II, where the neutron momentum distributions without FSI are identical to those in fig||. 

We observe that the momentum distributions are rather insensitive to the structure of the neutron-core subsystem 
although the distributions are narrower for a higher s-wave content. This is due to the fact that the final state 
interaction (neutron-neutron interaction) is identical in all the cases and therefore independent of the characteristics 
of the neutron-core interaction. Thus, the momentum distributions are sensitive to the 10 Li structure only through the 
n Li wave function. On the other hand, the effect produced by the final state interaction is significant, especially for 
a low- lying p^-resonance in 10 Li. Comparison with the experimental width for the neutron momentum distribution 
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in core break-up reactions again suggests a low-lying virtual s-state in 10 Li, see table II. 

C. Momentum Distributions for Finite Core Spin 

Up to now we have neglected the spin dependence of the neutron-core interaction or equivalently assumed that 
the spin of the core is zero. Any realistic calculation should of course assume that the spin of 9 Li is 3/2, and the 
neutron-core potential should include a spin dependence splitting the two s-states in 10 Li into states of total angular 
momenta 1 and 2. To do this, we have in the neutron-core potential included a term proportional to s„ • s c , where s n 
is the spin of the neutron and s c the spin of the core. For simplicity the spin splitting term has been introduced only 
in the s-wave, see eq.(|40|). 

The V s parameter has been chosen to be —7.14 MeV, the spin splitting parameter 7 S is taken as —0.13 for the solid 
curve and as 0.13 for the dashed curve. With these values the solid curves correspond to a situation where the lowest 
virtual s-state in 10 Li has spin 1 and an energy of 50 kcV, while the second s-state has spin 2 and energy 340 keV. For 
the dashed curves the lowest virtual s-state in 10 Li has spin 2 and energy 97 keV, and the s-state with spin 1 has an 
energy of 460 keV. In both cases, the average position of the virtual s-states is about 230 keV. In these calculations 
we have taken Vi=i = —19.14 MeV and V so — 6.97 MeV, placing the lowest p-resonance in 10 Li at 0.5 MeV. Final 
state interactions are included. 
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FIG. 5. One-dimensional neutron momentum distribution for a neutron removal process in n Li fragmentation. Solid curves: 
The virtual s-states in 10 Li are at 50 keV for s x — 1 and at 340 keV for s x = 2. Dashed curves: The virtual s-states in 10 Li 
are at 460 keV for s x = 1 and at 97 keV for s x = 2. For each case, the lower, the middle and the upper curves are obtained 
assuming that the spin of the final state system is 1, 2 and undefined, respectively. 

The momentum distributions can now be divided into two parts arising respectively from angular momentum 1 and 
2 of the final two-body state. This amounts to a separation of terms with j x = 1 and 2 in eq. (|33"|) . In principle they 
can be measured independently. In figj^ we plot results for the one-dimensional neutron momentum distribution for 
a neutron removal process. We have computed and plotted three curves for each case. The lowest curve (normalized 
to 3/8) is the one-dimensional neutron momentum distributions assuming that the final two-body system has spin 1. 
For the next curve (normalized to 5/8) the final two-body state is assumed to have spin 2. Finally the upper curve 
is the sum of the two others and therefore normalized to 1. For a final state of spin 1, the widths of the momentum 
distributions clearly increase with the energy of the virtual s-state. Compare the lowest solid and dashed curves 
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with FWHM of 26 and 37 MeV/c, respectively. This trend is also seen for a final state of spin 2. Compare here the 
intermediate solid and dashed curves with FWHM of 36 and 30 MeV/c, respectively. 

We now assume that the spin of the two-body system in the final state is not measured. Then the distributions 
for the individual spins should be averaged with the statistical weighting already included in the normalizations. 
This amounts to include the s x summation in eq. (p3|) . The two cases in fig.|| resulting in the upper the solid and 
dashed curves can hardly be distinguished (FWHM=32 MeV/c). This is due to the fact that the statistically weighted 
average energy of the two virtual s-states in 10 Li in both cases is about 230 keV. The Pauli principle dictates that 
both s-states independent of splitting are equally occupied in 11 Li. The influence of the p-resonance at 0.5 MeV is 
insignificant and the s-wave content is always about 91%. Therefore the narrow distribution arising from one state 
must then be compensated by the broader distribution from the other s-state. 

Thus, when the spin of the final two-body state is unknown the momentum distributions are only sensitive to the 
average position of the virtual s-states and the s-wave content in the initial wave function. Of course this conclusion 
also holds for core momentum distributions and core break-up reactions. If the momentum distributions could be 
measured for a given spin of the final two-body state, information about the lowest lying virtual s-state of 10 Li could 
be extracted from our analyses in a very direct way. 

IV. REALISTIC CALCULATIONS FOR n LI 

The calculations we made in the previous subsection, even though a very simple 11 Li structure was assumed, 
indicate that the unbound 10 Li nucleus has a low- lying virtual s-state (table II). A low lying p-resonance clearly 
overestimates both the width of the core momentum distribution and the neutron momentum distribution in a core 
break-up reaction. Recent analyses suggest that a virtual s-state is present at around 50 keV |14j together with a 
p-resonance at 0.5 MeV |30| . 

We first adjust the strength parameters in the neutron- 9 Li potential, where Vi=i and V so are chosen to place the 
lowest p-resonance at 0.5 MeV. The spin-orbit parameter V so is used to vary the p-wave content in the neutron- 9 Li 
subsystem of the total n Li wave function. This is achieved by moving the higher-lying p-resonances up or down. 
Then the values of the parameters V s and "f s are chosen to give a low-lying virtual s-state and at the same time 
reproduce the correct binding energy and rms radius of n Li. The lowest virtual s-state in 10 Li is assumed to have 
spin 2, but choosing spin 1 instead does not produce significant changes of the numerical results. The interesting 
physical parameters are then energy and size of Li, the p-wave content and the positions of the lowest s- and p-wave 
resonances in 10 Li. 

In fig|j| we show one-dimensional 9 Li momentum distributions from 11 Li fragmentation. The energy of the lowest 
virtual s-state is 50 keV. Different p-state contents in the neutron- 9 Li subsystem have been considered, 4% (solid 
curves), 18% (short-dashed curves), 26% (long-dashed curves), and 35% (dot-dashed curves). The thick and the thin 
curves are the calculations with and without final state interactions, respectively. The experimental data correspond 
to longitudinal 9 Li momentum distributions from a fragmentation process of a n Li projectile of energy 468 MeV/u 
and 648 MeV/u colliding with an Al target plj |. As expected and observed in fig|l], the effect of final state interactions 
is not very big, due to the large mass of the core. However, even these small effects seems to be necessary to reproduce 
the measured core momentum distributions, which in fact are rather insensitive to the structure of the neutron-core 
subsystem. All the curves including final state interactions match the experimental data equally well. 

Deviations between measured and computed curves increase with momentum. Two different sources for this dis- 
agreement are obvious. First the three-body model does not include core degrees of freedom. Therefore, if the 
neutron-core distance is smaller than the core-radius of about 3 fm, or equivalently the neutron-core relative momen- 
tum is larger than about 65 MeV/c, the model is not applicable. Secondly, only contributions from nuclear break-up 
reactions are included in the sudden approximation and other processes are assumed to be negligible. This is in 
general believed to be very well fulfilled assumptions, especially for the longitudinal core momentum distributions 
considered here. A convincing example is the fragmentation of 6 He into an a-particle and a neutron, see 0]. Our case 
of n Li is expected to be more sensitive to possible additional contributions than this very well reproduced transverse 
momentum distribution of the a-particle. The good agreement for the major part of the momentum range in the 
figj^ is an aposteriori confirmation of the validity of the model assumptions, i.e. the accuracy of the three-body wave 
function, the treatment of the reaction mechanism and the incorporation of final state interactions. 
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FIG. 6. One-dimensional core momentum distribution for a "Li fragmentation reaction. The lowest virtual s-state and 
p-resonance in 10 Li have an energy of 50 keV and 0.5 MeV, respectively. The thick (thin) curves are the calculations with 
(without) final state interactions. The p-wave content in the neutron- 9 Li subsystem is 4% (solid curve), 18% (short-dashed 
curve), 26% (long-dashed curve), and 35% (dot-dashed curve). Experimental data are taken from H. 
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FIG. 7. Same as in figjf] when the energy of the lowest virtual s-state in 10 Li is 200 keV. The p-wave content in the neutron- 9 Li 
subsystem is 7% (solid curve), 19% (short-dashed curve), 29% (long-dashed curve), and 33% (dot-dashed curve). 
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When the energy of the lowest virtual s-state increases the core momentum distributions become broader as seen in 
figj^, where the virtual s-state energy is 200 keV. The computed momentum distributions then overestimate the width 
of the experimental distributions, especially the ones with a higher p-wave content in the neutron- 9 Li subsystem. This 
conclusion is perhaps only convincingly reached after both the visual impression from the figure and consultation of 
the numerical widths in table II. We can then consider 200 keV as an upper limit to the energy of the lowest virtual 
s-state. On the other hand, a lower limit to this energy is not provided by the core momentum distribution, which 
remains unchanged for virtual s-state energies below 50 keV. 
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FIG. 8. Two-dimensional neutron core momentum distribution from a core break-up reaction. The lowest virtual s-state 
and p-resonance in 10 Li have an energy of 50 keV and 0.5 MeV, respectively. The interpretation of the curves is as in fig|| 
Experimental data are taken from jlq |. 

Let us focus now on neutron momentum distributions, where the effect of final state interactions is larger (figs.^, |^, 
and |]). Since for core break-up reactions the final state nucleon-nucleon interaction is well known, this kind of process 
is a good test for the model and the method of including final state interactions. In fig.^j we plot two-dimensional 
neutron momentum distributions for core break-up reactions for the same 10 Li structure as used in fig.|| (virtual s-state 
at 50 keV and p- resonance at 0.5 MeV). The meaning of the curves is also as described in figJq. The experimental 
data are obtained in a fragmentation reaction with a C target and a n Li projectile of 280 MeV/u fl5| . 

We observe that final state interactions now are crucial for reproducing the experimental distribution. The cal- 
culation with low p-wave content is too narrow, indicating that some p-state content is needed in the neutron-core 
subsystem. When we gradually increase the energy of the virtual s-state from 50 keV to 200 keV the momentum 
distributions remain esentially unchanged (the variations are of the size of the width of the curves). This insensitiv- 
ity to the neutron- 9 Li structure is not surprising, since the remaining particles (two neutrons) constitute a different 
subsystem. 

Finally, we compute two-dimensional neutron momentum distributions for a neutron removal reaction. We first 
consider the case of a low-lying virtual s-state in 10 Li with an energy of 50 keV. The two-dimensional neutron 
momentum distributions are plotted in fig|| The meaning of the curves is explained in fig.^j. The experimental data 
are taken from a fragmentation reaction with a C target and a n Li projectile of 280 MeV/u |14[ . We observe again 
that final state interactions are decisive in order to reproduce the experimental distribution. Since now the final 
state interaction depends on the details of the neutron-core potential, the computed distributions show a stronger 
dependence on the p-state content in the 10 Li subsystem. The momentum distribution with a low p-wave content 
underestimates the width of the distribution, while any neutron-core interaction producing more than 35% of p-wave 
content in 10 Li overestimates the width. The best agreement with the experimental distribution is found for a p-wave 
content ranging from 20 to 30%. 
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FIG. 9. Two-dimensional neutron core momentum distribution from a neutron removal reaction. The lowest virtual s-state 
and p-resonance in 10 Li have an e nerg y of 50 keV and 0.5 MeV, respectively. The interpretation of the curves is as in fig.[]. 
Experimental data are taken from UM. 




10 20 30 40 50 60 

p 1 (MeV/c) 

FIG. 10. Same as in fig.^|when the energy of the lowest virtual s-state in 10 Li is 10 keV. The p-wave content in the neutron- 9 Li 
subsystem is 4% (solid curve), 18% (short-dashed curve), 25% (long-dashed curve), and 33% (dot-dashed curve). 
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Let us now investigate how these neutron momentum distributions are modified by a variation of the energy of the 
lowest virtual s-state in the neutron-core subsystem. In fig.|l^ we plot the same distributions as in figj| for a low-lying 
virtual s-state at 10 keV. The momentum distributions are now much steeper at low values of the momentum and the 
intermediate region is not well fitted. None of the distributions are able to reproduce the measured curve in the whole 
range of momenta. If we increase the p-wave content in the neutron-core subsystem to fit the intermediate region, we 
obtain too high values at large momenta. The best compromise would differ considerably more both at intermediate 
and high momenta than the compromise in fig|^. Furthermore, the deviation at very small momenta also seems to be 
larger although only based on essentially one measured point with fairly large error bars. Thus, we conclude that the 
virtual s-state energy in 10 Li must be larger than 10 keV. 




FIG. 11. Same as in fig.|9| when the energy of the lowest virtual s-state in 10 Li is 100 keV. The p-wave content in the 
neutron- 9 Li subsystem is 4% (solid curve), 19% (short-dashed curve), 27% (long-dashed curve), and 35% (dot-dashed curve). 

For energies of the virtual s-state larger than 50 keV the experimental neutron distributions distributions are also 
reproduced fairly well. In fig.jll] we show the same distributions as in figs]^ and [l^, but for a virtual s-state energy at 
100 keV. The agreement is improved from fig.|l^ and remains almost unchanged from fig.^. Even for higher virtual 
s-state energies the nice agreement is maintained. This is due to the lack of accurate experimental points at low 
momenta, which causes the uncertainty in the position of the origin of the distribution. A larger number of points in 
that region would further constrain the energy of the virtual s-state. 

Comparing figsj|]ll| we see that the best agreement with the experimental neutron momentum distributions is 
obtained for a p-wave content between 20 and 30%. The reason is that the behaviour of the distribution for large 
momenta is determined by the p-wave structure of the neutron-core subsystem. This structure is essentially the same 
in all these figures, i.e. the lowest p-resonance is placed at 0.5 MeV and the higher p-resonance considerably above. 
Furthermore, the neutron momentum distributions arc rather insensitive to the position of the p-resonances. A very 
large variation of these energies are needed to change the shape of the distributions significantly at large momenta. 

V. SUMMARY AND CONCLUSIONS 

We have in this paper given a detailed formulation of a general method to describe fragmentation of high-energy 
three-body Borromean halo nuclei on light nuclear targets. In the previous papers in this series we established general 
properties of the structure of three-body halo nuclei. The most detailed information about these structures can be 
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obtained from measurements of momentum distributions of the fragments from reactions of halo nuclei with various 
targets. The inherent mixture of the initial halo structure and the reaction mechanism requires a description where 
both ingrediences are included. The present paper is devoted to construct a suitable method to study the properties 
of fragmentation products from halo nuclei. 

The initial wave function must be available with a good precision and three-body computations are necessary. The 
reactions are described in the sudden approximation, where one of the particles instantaneously is removed from the 
three-body system while the other two particles are left undisturbed. This approximation is very well adapted for 
the weakly bound halo nuclei moving with high energy relative to the targets. The remaining two particles of the 
projectile move roughly in the beam direction while they continue to interact after the collision. This final state 
interaction is essential especially when low-lying resonances are present in the two-body system. The model must 
include these effects. 

We use the recently developed method to solve the Faddeev equations in coordinate space. The initial three-body 
wave function is then accurately determined. We first briefly sketch this method, which is particularly well suited for 
an accurate treatment of the large-distance behavior. The low-energy scattering properties described by the scattering 
lengths are the decisive quantities determining the large-distance behavior of the three-body system. Therefore we 
parametrize the two-body interactions in terms of simple spin-dependent operators and gaussian radial shapes. We 
require that the relevant measured s- and p-wave scattering lengths as well as the singlet s-wave effective range are 
reproduced by the nucleon-nucleon interaction. The neutron-core interaction is adjusted to reproduce the size and 
binding energy of n Li and selected properties of 10 Li. The initial three-body wave function can now be accurately 
computed. 

The halo nucleus reacts with a light target nucleus and one of the three particles is removed instantaneously. The 
final-state two-body wave function is asymptotically approaching the plane wave limit, but at smaller distances is 
influenced by the interaction between the remaining particles. We include this modification by a distorted wave 
approximation. This final-state two-body interaction is identically the same as the interaction entering in the com- 
putation of the initial three-body wave function. We maintain this consistency throughout the paper. 

The model is now completely defined and we proceed to compute the transition matrix element which describes 
the probabilities for the various break-up processes. We integrate away analytically or numerically all non-observed 
quantities and end up with expressions for the observable momentum distributions of the remaining particles. Both 
one- and two-dimensional distributions are computed. To compare with the experimental results we tranform the 
particle momenta to the center of mass of the three-body halo nucleus. This is a straightforward but technically 
rather lengthy procedure. 

Numerical results are presented to illustrate general aspects. We first assume a spin-less core and vary the s and 
p-state properties of the neutron-core system. The physical parameters are the energy of the lowest virtual s-state, 
the energy of the lowest p-resonance and the p-wave admixture in the neutron-core subsystem of the three-body wave 
function. The size and binding energy of the three-body system is left unchanged at values corresponding to Li. 

Three different momentum distributions can be calculated, i.e. (i) the core momentum distribution after neutron 
removal, (ii) the neutron momentum distribution after neutron removal and (iii) the neutron momentum distribution 
in a core break-up reaction. They each carry different information. For (i) the effects of final state interactions 
increase and the widths of the distributions decrease with increasing s-wave content of the neutron-core subsystem. 
The centrifugal barrier is crucial for this behavior. For (ii) the final state interaction is crucial and the momentum 
distributions strongly depend on the structure of the continuum spectrum of the neutron-core subsystem. The lower 
the virtual state or resonances the larger the effect. For (iii) the well known neutron-neutron final state interaction 
has a significant effect on the distributions while the properties of the neutron-core subsystem are less important. 

The effects of finite core spin are necessary for realistic comparison with experimental data for halo nuclei. The 
resulting spin splitting of the lowest virtual s-states in the neutron-core system seems to produce quantitative and 
perhaps even qualitative differences. However, in most cases only the statistically weighted average energy is important 
due to the constraints from the Pauli principle, the parity conservation and the known energy of the three-body 
system. One exception occurs when the two-body final-state total angular momentum is measured in addition to the 
momentum distributions. Then the angular momentum of the relative state is determined and only one of the s-states 
contributes. The large difference, due to final state interactions, between the momentum distributions from the spin 
split states can in this way be observed. 

Finally, we perform realistic calculations for 11 Li and compare with the avalaible experimental data. Since the final 
state neutron-neutron interaction is well known, core break-up reactions constitute a good test of the method used to 
incorporate final state interactions in the description of the process. The uncertainties involved in the calculation are 
therefore reduced. The other type of process, neutron removal reactions involve the structure of the unbound system 
10 Li ( 9 Li +n), that simultaneously is crucial for the structure of the n Li three-body projectile. 

The comparison with the experimental data permits us to conclude that the final state interaction is an essential 
ingredient in the computation of the momentum distributions, even for core momentum distributions. The p-wave 
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content in the neutron- 9 Li subsystem has been concluded to be between 20% and 30%. From the longitudinal core 
momentum distribution and the radial neutron momentum distribution for a neutron removal process we conclude that 
the energy of the lowest virtual s-state in 10 Li must be between 30 keV and 200 keV. To improve the uncertainties 
in these conclusions, more accurate experimental data for neutron momentum distributions for low values of the 
momentum would be very helpful. 

In conclusion, a consistent model has been developed to analyze momentum distributions of the particles after 
high-energy fragmentation reactions of three-body halo systems. Final state interactions are essential. Application 
to 11 Li and comparison with available data led to rather severe constraints on the properties of both 11 Li and the 
unbound 10 Li nucleus. 
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APPENDIX A: HYPERSPHERICAL AND HYPERANGULAR COORDINATES 

We consider a system of three particles with masses rm and coordinates r 4 (i = 1, 2, 3). The Jacobi coordinates are 
defined as 

x fe = a l3 ( Ti -rj), y k = a {ij)k f > ' - r fc J , 

1/2 / , / xi/2 Al 



1 rrii rrij 
m mi-\-mj 



where {z, j, k} is a cyclic permutation of {1, 2, 3}, and m is an arbitrary normalization mass. 
From these definitions we find the momenta and \<L y k to be: 
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where p i7 p^, and p k are the momenta of the three particles. When these momenta are referred to the center of mass 
of the three-body system (p^ + p^ + p k = 0) we get 

p % = Ojjkxfc + - ' a {tj)k kvk , (A4) 



m 



rrij + m 

'3 

The hyperspherical variables {p, cti, Vt X i, Vt yi } are defined as 



Pj = -aijk xk + a {ij)k \t yk , (A5) 

3 

Pfc = ~(Pi + Vj) = -a(ij)kKk ■ (A6) 



2 i „,2 



yj , a t = arctan(xi/^) , (A7) 



where £l X i and Q y i define the directions of Xj and y^. 

Analogously we define the hyperspherical variables in momentum space {k, a Ki , flk x i^k v i} as 

K = \J k li + k li - a « = arctan(fc„/fc yi ) , (A8) 

where f2fc x i and Qui define the directions of k xi and k yi . 
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